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Using the iV-patch renormalization group method we investigate the flow of the quasiparticle 
weight in one-dimensional, weakly two-dimensional and fully two-dimensional Hubbard models. 
In one dimension we reproduce the Luttinger exponent that describes the disappearance of the 
quasiparticle peak towards lower scales. Further we analyze effects of the band structure, interchain 
hopping and the flow of the interactions in quasi-one-dimensional models. For the two-dimensional 
case we study how the suppression of the quasiparticle weight affects the flow to strong coupling. We 
find that the flow to strong coupling remains essentially unchanged. This strengthens the evidence 
for d-wave superconductivity in the weakly repulsive Hubbard model. 



I. INTRODUCTION 

Renormalization group methods have contributed in many ways to the understanding of the low energy properties 
of interacting fermion systems. Early work using these methods focused on one-dimensional systems 1 ; mathematical 
studies have been done for weakly coupled systems in one 2 and two dimensions 3-7 . General properties of interacting 
electrons were discussed in 8,9 . The exact RG schemes were used to prove the existence of the Landau- Fermi liquid in 
two dimensions under certain conditions 5-7 . 

In the course of the last few years, the method has been applied to highly anisotropic two-dimensional models 
like the Hubbard model on the square lattice and variants thereof 10-18 without any crude simplifications of the 
dispersion relation and phase space. The approximate RG schemes used in these studies are derived from exact 
renormalization group equations 5 ' 13,19-22 that can be applied to wide variety of models and that allow a clearer view 
on the approximations necessary in practical calculations. 

In these approximate studies in two dimensions, the flow of the self energy has up to now mainly been neglected, with 
the exception of the following studies. For the t, t' Hubbard model the Fermi surface flow was included dynamically 
in the RG flow 12 (see the Appendix of that paper). It was found that at a fixed density close to the van Hove 
density, t' gets reduced in the flow. The quasiparticle scattering rate was estimated in the t, t' Hubbard model 15 . The 
quasiparticle weight was calculated for the special case of the half- filled band by Zanchi 16 . He found a substantial 
reduction of the weight at the saddle points (0,±7r) and (±7r, 0). 

In the following, we discuss the question of selfenergy effects in a wider context, to bring out the points relevant 
for this paper. The selfenergy shifts the Fermi surface and it changes the Fermi velocity and the quasiparticle weight, 
as well as the curvature of the Fermi surface and the quasiparticle scattering rate. For weakly coupled systems with 
curved Fermi surfaces and without van Hove singularities, detailed phase space arguments were applied to prove that 
the curvature of the Fermi surface develops no singularity, the quasiparticle weight remains close to one and the 
correction to the Fermi velocity is small 4-6 . In particular, this flow does not change the instabilities of the Fermi 
liquid towards a superconducting state, which shows up as an initially marginal growth of the four-point function 
that leads to a divergence at a scale related to the critical temperature. Moreover, the old question of continuity 23 
of the results of perturbation theory in the limit T — > was solved: it was shown in Ref. 24 that the problem of 
. £h ' " anomalous diagrams" in Ref. 23 arose only because of an insufficient renormalization procedure (which amounts to 
the attempt of describing the shift of the Fermi surface merely by a shift in the chemical potential). 

The deformation of the Fermi surface was calculated by Halboth and Metzner in perturbation theory in the Hubbard 
model 25 ; they found that there is a density p ss 0.6, below which the Fermi surface tends to become rounder and 
above which it tends to become more squarelike. We note that this density value can be understood geometrically as 
the smallest one where Umklapp processes first contribute to the two-loop selfenergy; at this density the singularity 
analysis of Ref. 4 changes as well. Recently it was pointed out that a Fermi surface close to the van Hove points is in 
principle unstable towards a spontaneous breaking of the fourfold symmetry 26,27 . 

A curved Fermi surface away from the critical points of the dispersion relation is a necessary condition for Fermi 
liquid behavior. In presence of van Hove singularities and in situations with small curvature, such as very close to half- 
filling at t' = and in quasi-one-dimensional systems, the selfenergy can be expected to be relevant. The prototypical 
case of a perfectly flat Fermi surface is (apart from half-filling in two dimensions, where the situation is further 
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complicated by the presence of van Hove singularities) the one-dimensional fcrmion system. Simple descriptions of 
the flow towards the Luttinger liquid within the Wilsonian RG concept have been given for example by Metzner et 
al. 28 or the recent comprehensive review by Bourbonnais 29 . A mathematical proof that the flow of the quasiparticlc 
weight leads to the Luttinger exponents was given by Benfatto et al. 2 . Recently Busche et al. 30 rederived the flow 
of the quasiparticle weight Z in the ID Luttinger model and even obtained the spectral function of the Luttinger 
fixed point. Thus the functional RG schemes give a detailed description of Luttinger liquid physics. One important 
advantage of this method is that it neither requires exact scaling laws nor relies on the peculiarities of one spatial 
dimension, and that it applies to a wide class of model Hamiltonians. 

From the RG point of view, the existence of the Luttinger liquid fixed point hinges on the cancellations to all orders 
between the particle-hole terms at 2Uf and the particle-particle terms at 0. This cancellation is there already in 
lowest order, where the wave function renormalization still does not enter. It is a one-dimensional phenomenon: in 
higher dimensions, it is absent, and the interactions almost always flow to strong coupling, i.e. diverge at a positive 
energy scale. Already before this point, the RG flow in the fermionic variables breaks down. This signals the possible 
opening of a gap in the fermionic excitation spectrum, the simplest case being superconductivity. A key question 
is if the selfenergy effects, when fed back into the flow equations of the four-point interaction, can remove such 
instabilities altogether. In particular, a strong suppression of the quasiparticle weight might prevent a flow to strong 
coupling, and thus instabilities with ordering tendencies, at all. For curved Fermi surfaces, this has proved not to be 
the case 4-6 , but the situation is less simple for quasi-lD systems 31 ' 32 and models with van Hove singularities on the 
Fermi surface 16 ' 33 ' 34 . 

In the present article we address this question using the numerical iV-patch implementation of the RG scheme for 
the 1PI functions 13 . We incorporate the second order flow of the selfenergy in a way similar to the one-dimensional 
study by Busche et al. 30 , but we also include the flow of the renormalized interactions, which is relevant in more than 
one dimension. 

The main results are as follows. Our calculation correctly reproduces the suppression of the quasiparticle weight 
and the exponent of the Luttinger model in the one-dimensional Hubbard model. In quasi-lD systems with curved 
Fermi surfaces the flow to strong coupling remains almost unaffected by renormalization of the quasiparticle weight, 
in the parameter regime we study. The same holds for the fully two-dimensional system: the suppression of the 
quasiparticle weight does not prevent the flow to strong coupling. This implies that the (i-wave superconductivity 
found in the weakly repulsive Hubbard model 10-12 ' 14 close to half filling is not destroyed by the inclusion of these 
selfenergy effects. Another piece of information that is provided by our scheme is the anisotropy of the quasiparticle 
weight around the Fermi surface. At temperatures above the flow to strong coupling this information can be used 
for a Fermi liquid description of the normal state. Our finding is that in the simple Hubbard model the quasiparticle 
weight varies only mildly around the Fermi surface and no anomalous effects should be expected based on that. 

II. THE METHOD 

The renormalization group equations for the 1-particle irreducible (1PI) vertex function have been described in 
several works 12 ' 13 ' 21 ' 35 . The basic input for these equations is to assume that the quadratic part of the original action 
depends continuously on a certain parameter, here called A. The RG equations then describe the evolution of the 
vertex functions of the system when this parameter is varied. In most applications A is taken to be an infrared energy 
cutoff (or, similarly, momentum shells arc integrated out); under certain conditions however it may be appropriate to 
use the physical temperature 17 . Although the latter choice is intuitively appealing and preferable when investigating 
if ferromagnetic instabilities occur 17 , we will use here the conventional IR cutoff scheme, because it is simpler and 
requires less numerical effort in the iV-patch implementation, and because here we shall only study parameter regions 
where ferromagnetism was found to be absent 17 . 

In the IR cutoff RG the quadratic part Q of the action, in our case Q(iuj, k) = iu — e^, is supplemented with a 
cutoff function XA(k), 

T^2^2-ijj(iuj,k)Q(iuj,k)tl;(iuj,k) — ► T ^ ^ ^(iw, k)Q(iw, k)i/j(iw, k) Xa^W ■ 

i u k ^ k 

We take xa(&) = K(e^ — A) + K(—e? — A) where K(x) is a fixed smooth function that increases from to 1 in a 
small interval of length 2rj around x = 0. Later we shall take the limit 77 — 0, which gives a step function cutoff. The 
function K and the stepfunction limit are discussed in the Appendix. Thus Xa(^) 1 is unity for band energy |eg| > A 
but gets very large for |e^| < A, thus strongly suppressing low energies and effectively restricting the integration over 
the fermionic modes to those above A. Taking the derivative with respect to A one obtains a functional RG equation 
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for the generating functional of the connected correlation functions and its Legendre transformation, the effective 
action. An expansion in the fields yields an hierarchy of coupled differential equations for the 1PI m-point vertex 
functions. The hierarchy is infinite because the derivative of the m-point vertex depends on the (m + 2)-point vertex. 
We truncate the hierarchy by setting the 1PI 6-point function equal to zero, and are left with the two equations 
presented graphically in Fig. 1. The first equation determines the flow of the self energy Ea(«w, k) and the second 
that of the two-particle interaction, i.e. the four-point vertex. The barred internal line in Fig. 1 denotes a single scale 
propagator 

S A (iuj, k) = -G A (iuj, k)^-Q A (iuj, k)G A (iuj, k) = — — K^4\ XA ^) i W 

OA {iu _ e - - xa(«O s a(^, k)Y OA 

while the internal line without bar stands for the interacting one-particle Green function at the given scale, 

G A {iu, k) = =r . (2) 

iu - t£ - XA(k)T, A (uv, k) 

In this paper we consider the limit of a sharp cutoff, K(x) = 9(x). We discuss in the Appendix how this limit is 
taken; the result is 

GK{i»& = - dA{ t\. A (3) 

iu - ej: - Sa(iw, k) 

and 

S A (iw, k) = '- \ k — = 5- ^rO A (k) (4) 

iu) — eg — Y, A (iu), k) iu) — e^ — T, A (iijj,k) OA 

with A (k) — 0{e^ — A) + 0{—ej: — A). When inserting these expressions in the flow equations, undefined quantities 
like 6{E — A)S(E — A) can appear at special values of the shift momentum q (e.g. at q = 0). The sharp cutoff limit of 
all such terms is, however, well-defined and unique. We discuss in some detail in the appendix the limiting rule for 
avoiding these (seeming) ambiguities. 

By spin rotation symmetry, the coupling function is 

7i^ 2 s3S4( fc i' fc 2,fc3,fc4) = S Sl . S3 S S2 . Si V A (k 1 ,k2,k 3 ) - S S2tS3 S SuS4 V A (k 2 ,k 1 ,k 3 ) 

(with £)4 fixed by momentum conservation in terms of fci, k% and k 3 ). Here ki = (wj, ki) includes frequency and spatial 
momentum. 

In terms of an expansion in the bare coupling function, the above-mentioned truncation is exact up to second 
order. However, the RG provides more than just a second-order calculation: the evolution of the interaction and 
the selfenergy is continually fed back into the RG differential equation. This effectively sums up contributions from 
arbitrarily high order and thus produces a scale-dependent resummation of perturbation theory. We note that the 
above-mentioned truncation does not correspond to an expansion to a fixed loop order: the flow equations appear 
to be one-loop, but they also take into account two-loop effects by iteration. Since we are not doing a semiclassical 
expansion in orders of h, this is not a problem by itself. The relevant question is when higher orders significantly 
change the flow. They certainly do so when the coupling functions get too large; a detailed discussion of what " large" 
means is given in Ref. 13. 

By power-counting arguments, the leading part of the flow is given by the interaction processes between particles 
close to the Fermi surface, that is, with frequency close to zero and spatial momenta on the Fermi surface. Thus 
we drop the frequency dependence of the four-point vertex, so that V A (ki, fe, £3) = V A (k\, k 2 , ks), and calculate 
V A (ki, fe, £3) with wavevectors k\,ki and £3 on the Fermi surface. We discretize the remaining angular dependence 
by introducing patches around N wavevectors on the Fermi surface and approximate the coupling function by its 
0(N 3 ) values at the patch wavevectors. This device of angular patches (or "sectors") has been used in various works. 
It is the basic method for mathematical studies of 2D fermion systems, where N is chosen as A -1 / 2 at scale A, i.e. 
it increases at low energies, as one approaches the Fermi surface. The strong constraints implied by momentum and 
energy conservation in two dimensions then yield an intrinsic 1/N expansion 3 ' 8,9 . As a concrete calculational tool 
this A^-patch technique was first used by Zanchi and Schulz 10 and then in a number of subsequent works 11,12,14-18 . 
The projection to the Fermi surface suffices to obtain the leading coupling function flow but it would be incorrect to 
assume that the full IPI four-point function is well approximated by the projection to the Fermi surface. To obtain 
observables, one has to study the flow of response functions 10-12 , or, more generally, the flow of functions with general 
external momenta. 
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FIG. 1. Flow of two-point and four-point vertex. 
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FIG. 2. The coupling function Va(&i, fe, fe). The spin indices belonging to wavevectors k\ and £3 (and &2 and ki) are the 
same. 

The quasiparticlc residue that we want to calculate is determined by the frequency dependence of the selfenergy. 
Because we do not use a cutoff in the frequencies, we would in principle have to calculate the self energy over the 
whole frequency axis. However, the singularity of the propagator is at zero frequency, so the behavior of the selfenergy 
near to zero frequency is the most important. We calculate the frequency derivative at zero imaginary frequency and 
at the Fermi surface as 



cL£ A (0,£ F ) 



S A (i7rT,fe F )-S A (-^T,fc F ) 
2ttT 



(5) 



where the denominator is the difference between the two smallest Matsubara frequencies u>±i = ±nT (here kp denotes 
the projection of k to the Fermi surface). We then determine the quasiparticle weight as 



and write 



Z A (k) = l + id u E A (0,k F ) 



g a (iuj, k) = e A (k) = z A (k) c A (k) 



S A (iuj, k) = ^rO A (k) = Z A (k) ^rC A (k) . 



(6) 

(7) 

(8) 



The above involves a number of approximations. As already discussed, we have approximated the quasiparticle weight 
for a particle at band energy by the quasiparticle weight on the Fermi surface at the RG scale A. This effectively 
reduces the total weight to Z A because high frequency modes also get this weight. Since we are interested in the low 
energy properties, this plays no important role. We have also neglected the Fermi surface shift Xa(0, kp). Moreover, in 
writing (7), we have assumed that the correction to the fc-dependence of the dispersion relation (involving V£a(0, fc F )) 
is of the same order of magnitude as Z A so that we can effectively write in the denominator once we have put Z A 
in the numerator of the propagators. This is correct in the one-dimensional case because there, the dependence on 
fco and on k is similar; in 2D it is a further assumption. An investigation of this assumption is underway, as well as a 
more detailed calculation of the Fermi surface shift extending that of Ref. 12. The quasiparticle scattering rate in the 
2D case was considered in Ref. 15 and turned out to be reasonably small until the couplings exceed the perturbative 
range. 

As indicated above, the frequency dependence of the four-point vertex can be dropped in determining the leading 
part of the flow, but not from the full four-point function. In particular, it is obvious from the form of the equation 
for the selfenergy that the external momentum does not enter the internal line, hence the frequency dependence of the 
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selfenergy must come from a frequency dependence of the four-point function. Again we employ an approximation 
that is exact to second order - like in Ref. 15 we reinsert the integrated form of the equation for the four-point 
function in the selfenergy equation. This gives a two-loop diagram, in which we now again approximate the two 
appearing full four-point functions 7^ by coupling functions Va< ■ The right hand side of the differential equation 
for the selfenergy is now nonlocal in A: the change of four-point vertex and selfenergy at scale A involves four-point 
vertices and selfenergies at scales A' > A. In the next section we will compare this nonlocal equation with its local 
approximation where four-point vertex and self-energy at scale A' are approximated with the ones at A. 

This iteration of the solution of the RGDE is a general method that can also be applied to the full system of 
equations: setting the 1PI 2m-point function equal to zero means restricting the 2m-point function to tree diagrams 
in the bare interaction. If the bare interaction is a four-point interaction, these trees are of order r = to — 1 in the 
four-point coupling. The resulting truncated system of equations can be solved by iterative substitution as above. 
Thus the truncation of setting the 1PI 2m-point function equal to zero is exact to order to — 1 in the coupling, and 
indeed we will reproduce the Luttinger liquid exponents to second order in our calculation. By the same method one 
can get higher order corrections to the Luttinger exponents if one truncates at a higher m. 

We note that this method to reconstruct higher-loop contributions to the selfenergy is less approximate than the 
scheme used in Ref. 16. The latter assumes that the forward scattering occurring in the one-loop diagrams for the 
selfenergy can be decomposed in parts that only depend on the total momentum or the momentum transfer (which may 
be less accurate for anisotropic systems) and interchanges scale and frequency dependences of the scattering vertex in 
a non-generic way. The reconstruction of the two-loop contributions as described above renders these approximations 
unnecessary. 

There are two distinct one-loop diagrams for the flow of the self energy (the Hartree and the Fock diagram) and 
5 one-loop diagrams (one particle-particle and four particle-hole diagrams) for the coupling function (see Figures 
4 and 5 in Ref. 13). Therefore we get ten diagrams from inserting the flow of the coupling function into the two 
one-loop diagrams for Sa, shown in Fig. 3. Four of these diagrams are of one-loop x tadpole-type, hence frequency- 
independent in our approximation and thus do not contribute to Z A . Thus we are left with six two-loop diagrams for 
(d/dA)d u ,S A (0,k). 
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FIG. 3. Second order diagrams for the self energy. The diagrams correspond either to tadpole x one- loop contributions or 
real two-loop terms. We only keep the logarithmically divergent two-loop diagrams. One internal line is at scale A, indicated 
by the bar in Fig. 3, and one of the other two lines corresponds to a single-scale propagator at scale A' that is being integrated 
over. The third line is an interacting Green function with support on modes above A'. 



In condensed notation the RG equations for the selfenergy are 



= y W * Sa = dA' V A < * S A , * G A > * V A ^j * S A 



(9) 



Here the star stands for the momentum integration and Matsubara frequency summation according to the standard 
diagrammatic rules (for details, see Rcfs. 12, 13). Similarly, the equation for the coupling function reads 



dV A 
OA 



= V A * S A * G A * V A . 



(10) 



In the following we consider the renormalized vertex and coupling function 
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7i 4) SlS2S3 , 4 (fcl^2,fc 3 ,fc 4 ) = [^A(fcl)^A(fc 2 )^A(fc 3 )^A(fc4)] 1/2 7K lS2S3S4 (fcl,fe,fc 3 ,fc4) (11) 

V A (k u k 2 ,k 3 ) = [z A (fc 1 )Z A (fc 2 )Z A (fc 3 )Z A (fc 4 )] 1/2 V A {k u k 2 ,k 3 ) (12) 

They describe the interaction between the quasiparticles with renormalized amplitude Z A at the corresponding energy 
scale. If the limit A — > can be taken and Z A remains nonzero, the quasiparticle scattering amplitudes of the emergent 
Landau Fermi liquid can be read off eq. (12). 

The flow equations can now be expressed in terms of 7 and Z A and the Ca defined in (7) as 

^(k F ) = -idJ^*C A )(0,k F ) (13) 

(here the dot denotes Ad/dA) and 

4 Z\(k-) 

7(fci, fc 2 ,fc 3 ,fc 4 ) = (7 * Ca * Ca * 7)^1^2, fc 3 , fc 4 ) + -V 7(^1,^2, fc 3 , fc 4 ) (14) 

i=i 2 ^a(^) 

In an expansion in terms of the bare coupling function, the last term on the right hand side of (14) is of third order. In 
the Luttinger model, this term cancels exactly against the third order contribution to the coupling function flow 2 ' 30 . 
Because of the cancellation between the particle-hole and particle-particle terms in the first term, the coupling in 
the Luttinger model does not flow in this approximation, as it must be. In higher dimensions we also have to drop 
this term to remain consistent to second order. Keeping this term would require also taking into account the higher 
orders in the flow of the coupling function, that is, truncating only at 7^. 

Therefore, in the following we drop the second term in (14). As just discussed, this is consistent both in ID and 
in higher dimension and thus allows us to compare the flow in the Luttinger model with that in weakly and fully 
two-dimensional situations. An important feature of the equations is that the Z factors do not appear in the equation 
for the renormalized four point vertex any more, so that the flow of the renormalized four point function is the same as 
that of the (unrenormalizcd) four point function when the Z factor is dropped altogether. Therefore a small Z factor 
can only prevent a flow to strong coupling through (11). One conclusion of this work will be that - in contrast with 
the ID case - in the 2D system the suppression of the quasiparticle weight is weak as long as the coupling function 
remains in the perturbative range. The full two-loop flow of the coupling function is left for future work, but the 
value of Z obtained from our calculation will already give us a first estimate of the size of these two-loop terms. 

As discussed above, 7 is now replaced by the integrated equation (14), with V replacing 7 on the right hand side, 
in (13). We can now solve for Z in the one-dimensional case. If V is independent of scale and has initial value U, eq. 
(13) reads (in the local approximation described above) 

§ = U 2 a (15) 
with a the value of the two-loop contribution and U the coupling constant. This integrates to the power law 

A ^ 



z =(vJ (16) 

It follows that the two-point function has an anomalous decay exponent given by all 2 . This argument also shows that 
it is not necessary to include a rescaling of the momenta in the RG transformation to obtain anomalous dimensions, 
contrary to some claims 35 . 

In higher dimensions and for curved Fermi surfaces, the overlapping loop effect 4 implies that the two-loop integral 
is small, so that (13) reads 

Z _, tt2 fAlogA d = 2 



Z= W \A 5 d>3, 

the solution of which remains close to one for all A > if U is small. However, in that case, the coupling function 
is not constant because of the Cooper instability, so that eventually Z would also differ significantly from 1. On the 
other hand, to get a bound we only need to replace U by the maximal value of the A-dependcnt coupling function 
in (17). This shows that Z does not differ much from 1 before the couplings get of order 1/VA~, which is outside of 
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the range where the one-loop coupling flow is justified 13 . So Z plays no important role in the flow. Moreover, when 
restricting to temperatures above the critical temperature for superconductivity, the couplings remain small if they 
started out small, and the Z factor remains close to 1: the system is a weakly coupled Fermi liquid 5 ' 6 . 

When the Fermi surface contains van Hove points or flat parts, the Z factor can, however, get strongly suppressed. 
We investigate this in Sec. V. 



III. THE HUBBARD CHAIN 



In this section we apply our formalism to the one-dimensional Hubbard chain. We have already seen that when 
the couplings do not flow, the Z factor obeys a power law, so that our approximation does describe the Luttinger 
exponents. Thus the numerical calculations in this section serve mainly as a check for the numerical implementation 
but they also gives some insight on the magnitude of band structure effects. The tight binding dispersion along the 
chain is given by 



= — It cos k x 



(18) 



where the chemical potential ji determines the particle number. In the usual g-ology terminology 29 one considers 
four coupling constants g\,...,g$ with two incoming and one outgoing quasiparticle on the Fermi surface. For the 
Hubbard interaction U J2i n i,l n i,l & t the initial stage of the RG procedure they all have the same value gi = U. First 
let us consider the case away from half filling. In this case the one-loop flow of the interactions does not lead to 
strong coupling. The coupling g 3 involves scattering into states away from the Fermi surface and does not influence 
the low energy physics. <?4 does not couple to any logarithmically diverging channel and does not flow in the usual 
RG schemes. g\ is irrelevant and is driven to zero towards lower scales. (72 is somewhat reduced as the difference 
2.92 — gi is a constant in the one-loop flow. The simplest model completely neglects the flow of 52 and sets g± = 0. 
As discussed in the preceding section, and in accordance with Refs. 28 and 29, in the local approximation (described 
in the previous section) to the RG equations the quasiparticle weight is suppressed like 



(19) 



If we allow the coupling constants 171 and 52 to flow, this gets modified 29 to 



exp 





r A n 2 
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(20) 




A/t A/t 

FIG. 4. Flow of the quasiparticle weight in the Hubbard chain when the flow of the interactions is suppressed (left plot) 
or allowed (right plot). The thin solid line shows the analytical expectations described in the text. The thick line show the 
numerical results at zero temperature in the A- local approximation to the RG equations, while the thick solid lines are without 
the local approximation. The thick dashed-dotted line is for T = I0~ i t. All data are for fi = —t and U — 3t. 
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In Fig. 4 we compare these expressions with our numerical RG scheme. If we choose the same local approximation 
the agreement is excellent down to scales 10~ 5 t (we start the flow at Aq = t). We also observe some deviations in 
the initial flow due to the non-linearity of the dispersion relation. Below A w 10~ 2 t these deviations are negligible 
and could be absorbed in the powerlaws (19) and (20) by simple prefactors. The RG equations that are nonlocal in 
the scale yield lead to quantitative corrections. If we neglect the flow of the coupling function (left plot in Fig. 4), 
the non-local corrections slightly weaken the Luttinger exponent. In the case where the coupling function is allowed 
to flow (right plot), the non-local corrections suppress the Z factors further because the couplings gi and 52 are 
decreasing functions of the RG scale. Nonzero temperature T causes the expected saturation of the flow of Z at scales 
A < T. Note that in our case, where U = 3t, corresponding to g 2 / '(Aii' l v 2 F ) « 0.076, and the initial scale Ao = t, the 
suppression of Z is relatively weak. 



IV. EFFECTS OF INTERCHAIN HOPPING 



Next let us consider the case of nonzero hopping t± > between Hubbard chains. This leads to the dispersion 
relation 

€j[ = —2t cos k x — 2t± cos k y — fi . (21) 

The curvature of the Fermi surface (see Fig. 6) induced by t± prevents the cancellation of the one-loop flow of the 
interaction and cuts off the divergence of the self energy derivatives. Phase space considerations indicate that above 
this cutoff scale Aj_ set by the curvature the system will resemble a one-dimensional system, while below Aj_ it will 
behave like a two-dimensional Fermi liquid and maybe undergo a flow to strong coupling. Analytic expressions for 
the quasiparticle weight in coupled chains and the crossover from ID to 2D behavior can be found in the review by 
Bourbonnais 29 . They are based on a random phase approximation for the interchain motion and state that 

B 

^-^(b^width)^' ^ 

where the second factor describes the effective decrease of the interchain hopping by the renormalization of the 
quasiparticles on the chains. 9 is the exponent that describes the decay of Z\ on a single chain. In our case 8 -C 1. 
Thus we expect Aj_ w t±. 




10~ 4 1Cf 3 1(f 2 1Cf 1 
A/t 



FIG. 5. Left: Fermi surface for interchain hopping 0.3t and the 16 FS points used in the iV-patch calculation. Right: Flow 
of the most repulsive (solid line) and most attractive (dashed line) coupling constants. The repulsive scattering processes are 
largest between the almost nested flat sides of the Fermi surface (see arrow in left plot). The temperature is T = W~ 4 t, initial 
U = 3i. The minimal quasiparticle weight at A = T is 0.89. 

In absence of interchain hopping and Umklapp scattering, the one-loop flow of the coupling constants does not lead 
to strong coupling because the Cooper diagrams cancel the 2fcp-particle-hole diagrams that would otherwise lead to a 
nesting instability. This delicate balance gets destroyed by nonzero t± as particle-particle and particle-hole diagrams 
depend on the transverse wavevector in a different way. Partial nesting between almost parallel parts of the Fermi 
surface generates a growth of some coupling constants (see Fig. 5). Interactions between different layers of coupled 
chains can in principle induce spin density wave ordering with the corresponding wavevectors at finite temperature. 
With increasing t± « 0.5t and larger the Fermi surface becomes more and more rounded and sufficiently away from 
half filling the nesting features get washed out. 

The specific question we are interested in is whether the flow of the quasiparticle weight becomes significant in 
the scale range where these nesting features develop. This is important as many studies of potential instabilities in 
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quasi-lD materials neglect possible selfenergy effects. Here we only consider the change of the quasiparticle weight, 
but a thorough analysis would also include the renormalization of the band dispersion and lifetime effects. 

In Fig. 6 we compare the flow of the quasiparticle weight for three different values of the interchain hopping t±. 
The anisotropy of the Z- factor along the Fermi surface is very weak and can be neglected for these values of t' . Even 
for case of t± = the Z-factor is suppressed only weakly down to scales of 10~ 4 £ where the flow is cut off by the 
nonzero temperature 1 . For t± — 0.002i and t± — 0.005£ the suppression of Z is even weaker and gets cut off at scales 
~ tj_. Therefore, at least for the parameters considered here, the flow of the quasiparticle weight does not qualitatively 
alter the conclusions drawn from the one-loop flow of the coupling function. We expect that in most cases where 
the coupling function grows large at low scales the flow of Z can be neglected. In order to destroy the quasiparticle 
effectively one has to go to extremely low energy scales and temperatures. On the other hand, the suppression of Z 
towards lower energies may still be visible in tunneling or optical measurements. 
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FIG. 6. Flow of the quasiparticle weight in the Hubbard chain coupled by a interchain hoppings t± — (dotted line), 0.002i 
(solid lines) and 0.004£ (dashed line). The anisotropy of the quasiparticle weight over the Fermi surface is below the numerical 
precision. The temperature is T = W~ 4 t, initial U = 3t. 

V. THE QUASIPARTICLE WEIGHT IN TWO DIMENSIONS 

Next we study the 2D Hubbard model close to half filling. The one-loop flow of the coupling function has been 
analyzed extensively for zero 10 ' 18 and nonzero values 11,12 of the next-nearest neighbor hopping parameter t'. As 
discussed earlier, in our approximation the flow of the renormalized coupling function is unmodified compared to the 
earlier studies neglecting the flow of the quasiparticle weight. Here we focus on the case t' = and on the flow of the 
quasiparticle weight. We want to analyze a) how strong the suppression of the quasiparticle weight is and whether 
these two-loop effects get comparable to the one-loop terms that drive the coupling function; b) whether the Z-factor 
develops a distinctive anisotropy around the Fermi surface, e.g. at higher temperature, when the flow of the coupling 
function remains finite. 

First let us describe the TV-patch results shown in Fig. 7 for half filling and in Fig. 8 for electron density (n) = 0.9 
per site. We observe that the quasiparticle weight gets increasingly suppressed toward low scales where the coupling 
functions get large. The suppression is strongest close to the van Hove points. Thus the flow is consistent with a 
vanishing quasiparticle weight at low energies at least in the saddle point regions. Nevertheless we note that even if 
we follow the growth of the renormalized coupling function far out of the perturbative range, the suppression of Z 
remains rather weak. This shows that the two-loop effects diminishing the quasiparticle weight are small compared 
to the one-loop contributions driving the coupling function and that the one-loop approximation 10-12 for the latter 
may be justified in this scale range. 

A similar picture is found away from half filling. In Fig. 8 we show data for (n) = 0.87 and initial U = 2.25i. At 
low temperatures the flow goes towards a d-wave superconducting instability. The quasiparticle weight gets somewhat 
suppressed towards the instability, but the renormalized coupling constants flow off to strong coupling. Although the 
flow of the Z-factors is driven by the interactions, it still continues an approximately logarithmic decrease while the 



1 We refrain from going to lower scales as the numerics for many patches becomes very slow below these scales. 
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interactions seem to diverge with some power-law. This indicates that at these scales the two-loop contributions are 
much weaker than the one-loop terms that drive the flow to strong coupling. At higher temperatures, T > 0.2i the 
coupling function does not diverge, and we can follow the flow to lowest scales. The flow of the quasiparticle weight 
gets cut off at scales A ss T and remains close to unity over the full Fermi surface. Like in the quasi-lD cases studied 
in the last section the suppression of Z is only mildly anisotropic. 

Next we interpret these findings and put them into context with other works in the literature. The flow of the 
quasiparticle weight close to the saddle points was studied by Dzyaloshinski 33 for t' ^ in a system restricted to the 
saddle point regions, and more recently by Zanchi 16 for the half-filled case using a A-patch formalism related to ours. 
Dzyaloshinski 33 neglected the logarithmically divergent particle-hole diagrams with respect to the logarithm-squared 
particle-particle diagrams. Consequently in the repulsive Hubbard model his flow does not lead to strong coupling. 
Nonetheless the quasiparticle weight approaches zero at low scales, similar to the Luttinger liquid case. However, at 
least for our choice of parameters, the particle- hole diagrams always drive the flow to strong coupling, even if the 
nesting is weaker. Zanchi 16 argued that for the half-filled case and t' = the quasiparticle weight approaches zero at 
the saddle points at scales above the divergence scales of any susceptibilities, indicative of an anomalous normal state 
without long range order. The suppression of the quasiparticle weight is in qualitative agreement with our results. 
However we note that in the perturbative range, where the approach is valid, the suppression of Z is only weak. 
Therefore in our opinion no strong conclusions regarding the order of closely competing poles in the flow of Z and 
various susceptibilities can be drawn. In view of the weakness of the two-loop effects it appears to be more likely that 
for the half-filled perfectly nested case the barely attenuated growth of the coupling function leads to spin-density 
wave (SDW) order in the ground state where quasiparticles exist only above the SDW energy gap or at temperatures 
above the instability. This is corroborated by the fact that at temperatures above the runaway flow of the interactions 
(see thin lines in Fig. 7) the flow of the quasiparticle weight saturates at low scales and remains nonzero. Thus, 
although we cannot rule out the possibility, we do not find any evidence for a non-Fermi liquid or Luttinger-liquid 
like state at temperatures above the flow to strong coupling. 

A similar expectation regarding the one-particle spectrum holds for densities less than half filling: the growth of 
the renormalized coupling function out of the perturbative range indicates the opening of a gap at least on parts of 
the Fermi surface. In the case (n) = 0.9 per site we expect superconducting long range order of <i x 2_ J/ 2-symmetry. 
In our approximation the flow of the interactions is unchanged by the flow of the quasiparticle weight and develops 
a strongly dominant diverging pair scattering in the d x i_ y i -channel as discussed by several authors 10-12,14,18 . The 
weak suppression of the quasiparticle weight down to scales where the pair scattering is already much larger than the 
bandwidth shows the irrelevance of two-loop effects. Thus our results affirm the evidence for d-wave superconductivity 
in the Hubbard model at weak to moderate U. 

Note however that there are examples of one-dimensional systems where a one-loop flow to strong coupling does 
not lead to a (quasi-)long-range ordered ground state. In the half filled two-leg Hubbard ladder 36 the single particle 
spectrum is fully gapped, yet all spin- and charge correlation functions remain exponentially decaying. We pointed 
out 12 that in the 2D Hubbard model for nonzero t' the flow to strong coupling strongly resembles the one in the 
two-leg ladder and that this may indicate a novel strong coupling state for a certain parameter range of the 2D 
Hubbard model different from the t' = case studied above. 
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FIG. 7. Left: Flow of maximal and minimal coupling constants for the almost half-filled t' — Fermi surface (see inset) with 
chemical potential [i — — O.Olt, initial U = t and N = 32 points, at T = O.OOlf (thick lines) and T = 0.03t (thin lines). Right: 
Flow of the quasiparticle weight close to the saddle point (solid line) and in the BZ diagonal (dashed line) at T = O.OOlt (thick 
lines) and T = 0.05t (thin lines). The inset shows the small variations of Z around the FS for T — 0.0014 at the scale where 
the largest coupling constant is Wt (upper line), 20t and 30t (bottom line), i.e. already far out of the perturbative range. 
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FIG. 8. Left: Flow of maximal and minimal coupling constants for the t' = Fermi surface with electron density 0.87 per 
site, initial U = 2.25t and N = 32 points at T = O.OOli (thick lines) and T = 0.2t (thin lines). Right: Flow of the quasiparticle 
weight close to the saddle points (solid line) and in the BZ diagonal (dashed line). Again thin lines are for T = 0.2t and thick 
lines for T = O.OOli. The inset shows the 32 point on the Fermi surface used in the A-patch calculation. 



VI. CONCLUSIONS AND OUTLOOK 



The iV-patch renormalization group scheme enables us to trace the change in behavior of weakly correlated electron 
systems as a function of dimension, by interpolating the ratio of the hopping parameters in the two directions from 
zero to one. It reproduces the Luttinger liquid behavior in the one-dimensional models away from commensurate band 
filling on one side and the Fermi liquid state with its various infrared instabilities on the other side. The present RG 
calculation still involves a cascade of approximations but our results match well with established knowledge acquired 
by other methods. In one spatial dimension the quantitative precision of the RG scheme can be checked directly 
by exact numerical techniques: in the study of boundaries of Luttinger liquids the comparison with density matrix 
renormalization group calculations turned out to rather promising 37 . Biermann et al. 38 used extended dynamical 
mean-field theory (DMFT) to calculate the quasiparticle weight in coupled chains. Their findings for nonzero interchain 
hopping are qualitatively similar to the RG results, e.g. the suppression and also the anisotropy of the Z-factors are 
relatively weak along the Fermi surface (note however that the RG calculation for Z applies only at weak couplings, 
and higher order corrections might further decrease the Z-factors) . We believe that the true strength of the A^-patch 
RG lies in its flexibility and transparency. Due to these properties the method may open the way to an unified 
qualitative picture of low dimensional weakly coupled electron systems. 

A direct conclusion of this work is the partial justification of the neglect of the wave function renormalization in 
several previous renormalization group studies of the 2D Hubbard model 10-12,14 . At least for the parameters used in 
these studies and over the scale range where the perturbative scheme is expected to be valid, the quasiparticle weight 
remains nonzero and its flow does not affect the one-loop flow to strong coupling in the second order approximation. 
Similar statements can be made about the one- loop Fermi surface shift 12 and the quasiparticle scattering rate 15 . 
The physical picture that is suggested by these findings for the Hubbard model close to half filling is that the most 
plausible interpretation of the runaway flow of the coupling function is the formation of an energy gap at least on 
parts of the Fermi surface and not a gapless state resembling a Luttinger liquid. Due to the limited validity range 
of the weak coupling approach our results can not support a conjectured pseudogap 16 at temperatures above the 
antiferromagnetic instability in the half-filled Hubbard model. On the other hand our calculations show that the 
antiferromagnetic instability and also the d-wave superconducting instability found in earlier studies 10-12,14 persist 
if the flow of the quasiparticle weight is included. In particular the possibility of d-wave superconductivity in the 
Hubbard model is still debated 39-41 , however mostly at larger U. Our analysis underpins the positive answer for the 
weak coupling sector of the model. 

An analysis of the renormalization of the dispersion relation and of the higher-order flow of the coupling function 
will provide further checks for the existing results. Regarding further applications of the method we mention the wide 
class of quasi-lD conductors like Bechgaard or Fabre salts 42 . These systems show a variety of magnetically ordered 
or charge-ordered states and many undergo transitions into supposedly unconventional superconducting states. The 
hopping parameters of these materials are known to a good extent 43 and a weak coupling description is expected to 
make sense. Thus our method should be well applicable and useful to understand the main features of the phase 
diagrams of these materials and the interplay between one- and higher-dimensional physics. Some steps in this 
direction have been undertaken by Duprat and Bourbonnais 44 . 
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We have used an approximation scheme where all relevant contributions up to second order in the interactions are 
taken into account. This correctly reproduces the Luttinger liquid case, where the flow of the renormalized interaction 
vanishes also in second order 2 ' 30 . In those cases where the second order flow leads to strong coupling, higher order 
contributions will become large at some scale. The physical content of the strong coupling fixed points in two-loop 
g-ology studies 1 ' 29 is, however, rather questionable, as the fixed point coupling constants are typically larger than the 
bandwidth and thus outside the perturbative range. Moreover the information obtained by careful one-loop studies 
of ID systems agrees well with the results of numerical techniques or bosonization. Of course, our RG method can in 
principle be used later on to investigate the higher order flow of the interactions as well. For the time being, note that 
our present scheme at least provides estimates when higher order effects such as the flow of the quasiparticle weight 
become important. 

In addition we note that the Af-patch approach in its present form can be used to extract a precise Landau-Fermi 
liquid description of two-dimensional systems like the Hubbard model at temperatures above the flow to strong 
coupling, extending existing studies 45 ' 46 that do not consider the renormalization of the quasiparticle weight. In view 
of the somewhat harmless flow of the Z-factors in the two-dimensional models discussed above we do, however, not 
expect drastic changes of the results. 

ACKNOWLEDGEMENTS 

We thank Patrick Lee, Walter Metzner, Maurice Rice and Ashvin Vishwanath for valuable discussions. This research 
was supported in part by the National Science Foundation under Grant No. PHY99-07949. C.H. acknowledges 
financial support by the DFG (Deutsche Forschungsgemeinschaft). 

APPENDIX A: SHARP CUTOFF LIMIT 

We consider a general class of cutoff functions given as 

XA(p)=K(e p -A) + K(-e p -A) (23) 

where K is a smooth function that increases from to 1 in a small neighborhood of size rj > of zero. K is 
constructed as follows. Let h be a smooth function of x that vanishes for \x\ > 1 with the properties h(x) > for all 
x and h(x)dx = 1. Define H(x) = h(u)du and let K{x) = H{x/nf). As rj — > 0, K approaches the Heaviside 
function: lim^o K(x) = &(x), and thus becomes a "sharp" cutoff function. We note that we have to take the limit 
r\ — ► or else use this as a cutoff function only for A > rj because for A < rj, K no longer vanishes near zero. The 
following argument for taking the limit rj — > applies with trivial modifications also to the often taken choice of a 
Xa(E) that depends only on the ratio E/A. 

In the following we show that in the sharp cutoff limit, the propagators appearing in the RG flow equation are 
given by (3) and (4). Moreover, since the single scale propagator contains 5(±ep— A) and the full propagator contains 
9{±tp— A), formal expressions like 5{x)0{x) appear and we determine what they really mean in the limit. As it turns 
out, the correct rule is that one has to replace the cutoff function by a variable and integrate it from to 1. We 
discuss the integral 

I±(q) = J dpG A (q±p)S A (p) $ A (p, g ) (24) 

which is prototypical for the integrals appearing in the flow equation. The other cases appearing are that where G A is 
replaced by 1 (in the equation for the selfenergy) and when Ga is replaced by a product of Ga's (such terms appear 
in the flow equations for the higher m-point functions. The function $a is a product of vertex functions. 

Let us first consider the case without van Hove singularity. In this case, we can introduce a radial coordinate E = 
and an angular coordinate 9, so that p~p(E, 9) is now a function of E and 9. We also write p(E, 9) for (po,p(E, 9)). 
We have I± {q) =TJ2 Po J d9 (X+ + A_ ) with 

X+= dE J(E, 6) S A (p)G A (q ± P ) Hp, q) (25) 
Jo 

and A_ the integral over negative E. Because we restricted to positive E, Xa(p) simply becomes K(E — A). Inserting 
the definition of K in terms of H and changing variable from E to u = (E — A)/i], we get 
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X + = r duJ(A + u V ,9)- ^\ GA{q±p) <s> A{pjq) (26) 

J-A/n {Wo -A-U7]- H(u)E A (p(A + U7], d)y 

Here J is the Jacobian of the change of variables; at fixed E, its integral over 9 gives the density of states at E, 
and we inserted a 9 function to restrict the integration tou> — A/t]. Because we are at a positive temperature, the 
integrand is bounded by an integrable function. Moreover its limit as 77 — > exists pointwise. Thus, by dominated 
convergence, the limit rj — > exists and can be taken inside the integral, so that 

X+ = JZ dU J(A ' 9) (ipo-A^g^M))^ GA( * ± P) ^ q) (2?) 

Note that the limit of G\ is not a continuous function, but this does not matter for the exchange of limits. The 
function $a also has a limit as rj — > 0. This follows by the argument we just outlined and by iteratively solving the 
differential equation. Let q ^ 0. We can change variables to H and are left with 

X + = J { A,e)G A{q ± p{ A,e))^ M [ *H ^_~^^ ew (28) 

which corresponds to the statement made above that the cutoff function has to be replaced by a variable H that gets 
integrated from to 1. The integral gives (ipo — A — Xa(A, 0)) _1 so the final integration over 9 and p corresponds 

to a p-integral with a propagator S A + \p) = 5{e,p— A) (ip a — e^— Sa(p)) -1 . The contribution from X_ is, of course, 
similar, with S(ep+ A) replacing 5(ep — A). 

For the special case q = the energies and cutoff functions in Ga are identical to those in 5a, so that 

X + = J(A,6)$(p,q) fdH- -(*Po-A)ff 

Jo (m - A - H?: A (p ,p(A,6))) 3 

The integral over H now gives — \{ipo — A — £a)~ 2 . Note the extra factor 1/2. Let us compare this with the result 
of taking q — + in the above results. If we integrate over 9, set qo = and take the limit q — > 0, we also get only 
an integral over "half of the Fermi surface: the intersection of the level set La = {p ■ £p = A} (along which we 
integrate) with the support of ^(e p + g — A) (where Ga 7^ 0) is, in the limit q — > with q = q/\q\ fixed, equal to 

L A ~^ = {p e La : p • q < 0}, which is the half of the Fermi surface. 

Thus for the one-loop integrals, setting q = gives the same result as the average over q of the limits q — > with 
q fixed. The above argument also shows that seeming ambiguities in the limit are resolved by properly doing the 
LL-integrals. Looking at the details of the H integration is also necessary for q ^ if q happens to be a perfect nesting 
vector for the Fermi surface. 

If there is a van Hove singularity on the level set La, the Jacobian J(A, 6*) has a nonintegrable singularity at each 
point p* on La where Ve^* = 0. We cut out small neighborhoods of all such singularities by a smooth partition of 
unity. The above argument applies to the contributions from outside these neighborhoods. In a small neighborhood 
of each saddle point, we can now change variables from p to x, y so that — A = xy. The Jacobian j of this change 
of variables is smooth if the second derivative of e at p* is nondegenerate. The X + defined above now changes to 

X+= I dxdyj(x,y) V 'J 1 } 71 - x (30) 

J x 2 +y 2< r 2 [tp -A-xy-i: A (p ,p(x,y))H(r] L xy) z 

One can now scale in any way one wants; the resulting integral diverges logarithmically in rj. Thus at these A, the 
right hand side of the flow equation diverges, leading to an infinite slope in the solution as a function of the scale A. 
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